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^ ■ Abstract 

t^- \ The jump process introduced by J. S. Bell in 1986, for denning a quantum 

field theory without observers, presupposes that space is discrete whereas time is 

£NJ ■ continuous. In this letter, our interest is to find an analogous process in discrete 

time. We argue that a genuine analog does not exist, but provide examples of 

\q . processes in discrete time that could be used as a replacement. 
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One of the central challenges for "hidden variable" approaches to quantum mechan- 
ics, such as the de Broglie-Bohm pilot wave theory, is to provide an adequate account 
^— » ■ of relativistic quantum field theory. To address this, Bell introduced a jump process on 

a discrete lattice jU EH UJ] , intended to reproduce the quantum mechanical predic- 
tions for fermion number density in space. The same method can be used to generate 
stochastic trajectories for any discrete observable, both in field theory and in nonrel- 
ativistic quantum mechanics. For a discretized position observable in nonrelativistic 
quantum mechanics, Bell's process reduces to the de Broglie-Bohm pilot wave theory in 
the continuum limit [Tol ll4|. so it is a natural analog of this theory for discrete "beables" 
0- 

Although the "beables" in Bell's process are discrete, it still contains a continuous 
time parameter. However, there are several reasons for developing a discrete-time version 
of the process. Firstly, some approaches to quantum gravity are based on fundamentally 
discrete space-time structures, so a realist account of these theories along Bohmian 
lines would have to be fully discrete. Secondly, "hidden variable" theories, no matter 
whether they are realized in nature or not, can be useful for numerical simulations [T^ ITj. 
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visualizations [3113, bookkeeping |S., and obtaining better intuitions about quantum 
phenomena. Numerical simulations are discrete by nature, and a fully discrete theory 
may also be useful when dealing with quantum phenomena usually described in a discrete 
setting, such as those considered in quantum information and computation. Thirdly, 
Valentini has recently proposed that matter in quantum nonequilibrium, i.e. beables 
with distributions other than |\I/| 2 , if existant, may provide astonishing computational 
resources, enabling us to solve NP-complete problems in polynomial time. However, 
since classical analog computers can also outperform Turing machines if the continuous 
variables can be manipulated with perfect accuracy, this claim would be simpler to verify 
in a fully discrete model. 

In this letter, we highlight the difficulties inherent in discretizing Bell's jump process, 
and propose two concrete discretized processes that circumvent them and converge to 
Bell's process as the time step r goes to zero. Other possibilities exist, along the lines 
of recent proposals by Aaronson pQ, and these will be developed in future work. 

Bell's process is a Markovian pure jump process (Qt)tm on a lattice Q with rate for 
the jump q' — > q given by 

[llm{%\P{q)HPffl\%)] + 
<rt{q\q) = in - 



(* t |P(gO|* t ) 

where x + = max(x, 0) denotes the positive part of x 6 K, is the state vector of a 
quantum (field) theory, evolving in some Hilbert space J$? according to 

ih^± = H% , (2) 
dt 

H is the Hamiltonian, and P(q) is the projection to the subspace J^ q C Jif, where the 
J{? q form an orthogonal decomposition, = © geQ Jf? q - Relevant properties of Bell's 
process are that at every time t, the distribution of Q t is the quantum distribution 

(%\P(q)\%), (3) 

and that its net probability current between q' and q, a t (q\q')T'(Q t = q') — cr t (q'\q)F(Q t = 
q) where P denotes "probability," agrees with the quantum expression for the probability 
current, 

|lm (*\P(q)HP(q')\V) . (4) 

Since many constructions are easier in discrete time than in continuous time, one 
might have expected that there is an analogous Markov chain (Qt)terZ on Q with discrete 
time step r such that the probability ¥ t (q' — > q) for the transition q' — > q, i.e., the 
conditional probability F(Q t+T = q\Q t = q'), is given by a formula similar to (JTJ), with 
H replaced by a simple function of the unitary U defining the time evolution 

^ t+T = U% , (5) 
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and that one could arrive at this formula by a reasoning similar to the one leading to 
(HD from © and (gj), as given in Sec. 2.5]. 

However, this is not possible in any obvious way. The obstacle is that in the time- 
discrete case there is no obvious formula for the net probability current J(q, q') between 
q' and q, replacing (j3J) of the continuous case. Given an expression for J(q, q') in terms 
of P, and U, we could set 

which would define a Markov chain (QtjterZ whose probability current 

P t (g' -> q) P(Q t = q') - P t (g - q') P(Q t = q) (7) 

coincides with J(q, q') and whose distribution at any time t coincides with the quantum 
distribution (J5J), provided J(q,q') has the following properties: 

J(q, q') E R (8a) 
J(q',q) = -J{q,q') (8b) 

^j( M r<($iP( 9 ')i$) (8c) 

]T J{q, q') = (*\U*P{q)U\V) - (*|P(g)|*> • (8d) 

q'dQ 

Currents of the form ( |7| ). with transition probabilities © and distribution ©, have 
these properties by construction. Property (f5c|) expresses that no greater amount of 
probability can get transported away from q' than present at q', and (j8d|) guarantees the 
quantum distribution (jSJ) at the next time step. The obvious way of guessing a formula 
for J(q, q') is to start from one of the expressions 

(V\U*P(q)UP(q')\V) (9a) 
(*\P(q)UP(q')\*) , (9b) 

to multiply it by any numerical constant, to take the real or imaginary parts to ensure 
(JHaJ), and to anti-symmetrize in q and q' to ensure (|8b[). However, all expressions thus 
obtained generically violate (jSHJ), except for the anti-symmetrization of 2Re(l9*aTl. 



J{q, q') = \m (U*P{q)UP{q') + P(q')U* P(q)U- 

U*P(q')UP(q) - P(q)U*P(q')U) |^) , (10) 

which can violate (|5cjl (numerically we found 46 examples of such violations among one 
thousand randomly chosen U and ip in J$? = C 3 with fixed one-dimensional projections 
P(q) and P{q')). 
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However, a different reasoning leads to a process in discrete time that has some 
features in common with Bell's process. Choose H such that 



U = e~ lTH , (11) 

so that the evolution (J2J) generated by H is a continuation of the evolution (J3J) generated 
by U. (The degree of non- uniqueness of this choice is discussed later.) Then, consider 
Bell's process (Qt)tm in continuous time for this H. By restriction to just the integer 
times, we obtain a Markov process Qt '■= Qt for t G rZ. 

The process (Qt)terZ has the quantum distribution ® at every time. It is important 
for this that the two evolution laws (J2J) and (j5J) for \l/ lead to the same ^ t at every t that is 
an integer multiple of r. It makes no sense to ask whether the probability current of this 
process, F(Q t + T — q,Qt — q') — ^(Qt+r = q', Qt = q), agrees with the one prescribed by 
quantum theory, since, as discussed above, quantum theory does not prescribe a unique 
current in the discrete-time case. Note that in the limit r — > the process approaches 
Bell's process. This fact and the simple and straightforward construction of (QtjterZ 
suggest that this may be the closest one can get to an analog of Bell's process in the 
time-discrete case. 

The transition probability F t (q' — > q) = F(Q t + T = q\Qt = q') does not, however, 
possess a simple formula in terms of ty t , U, and P(-) analogous to ((H), only the following 
one: 

to+r to+r to+r 

Pto(V ->9) = 5Z Yl 5 i'm 5 4,in j / dt 2 --- dt n x 

n=0 q ,...,q n eQ tQ h tn _ t 

( T \ n 

x exp - / a s (Q\ } )ds Y[(?t k (qk\qk-i) , 

\ to ' k=1 

with cr s (q\r) given by ((TJ) and cr s (Q\r) := J2 q eQ a s(q\ r )- Eq. (fT2|) is a fact about any 
jump process in continuous time with jump rates o (applied here to Bell's process Qt)- 1 
The process Q is not completely determined by \l/o, U, and P(-) since H is not 
completely determined by ([lip, even though in many cases there may be a natural choice 

1 To get a grasp of Ijl2|l . begin with noting that <7 s (Q\r) is the total jump rate at time s in the 
configuration r. The probability that no jump takes place before time t, if the process starts at to 
in go, is exp(— J t cr s (Q\qo) ds). Thus, the probability that the first jump takes place between time t 

and t + dt is exp(— J {q <y s {Q\oa) ds) Ct(Q|<Zo) dt. The probability that the destination of the first jump 
is qi, given that the jump takes place at time t, is <J t (qi\qo)/a t (Q\q ). Conditional on that the first 
jump occurs at t and leads to qi, the distribution of the times and destinations of the further jumps is 
the same as for a process starting at time t in q t . Thus, the probability of a path q , . . . , q n with the 
fc-th jump between tk and tk + dtk and no further jump before to + r is the integrand of (|12|l times 
dt\ ■ ■ ■ dt n . Now add (respectively integrate) the probabilities of all ways the process can move from q' 
to q in the time interval [to, to + r], namely by means of n jumps at times t\, . . . ,t n with destinations 
qi, . . . ,q n . For a more detailed discussion of such probability formulas, see |10j and 



4 



of H. For example, if U has an eigenvalue e , then H may have as the corresponding 
eigenvalue any of the numbers - + —k with k G Z. More generally, for any self-adjoint 
operator S with spectrum contained in ^Z and commuting with H (in the sense of 
commuting spectral projections), H + S is another solution of (fTTj) for given U. A 
unique H could be selected by the additional condition that the spectrum of H be 
contained in (— -, -]. 

In the particularly simple situation \Q\ = 2, there does exist a time-discrete analog 
(QtjterZ to Bell's process. In this case, the expression (fTU|) satisfies (JHJ) and thus defines 
a process; in fact, the net probability current between the two configurations q' and q 
is already determined by the distribution (jSJ) and must be 

(%\U*P(q)U\%) - (V t \P(q)\%) , (13) 

since any increase or decrease can occur only by transitions from or to the other config- 
uration. Just as Bell's process has the smallest jump rates compatible with the current 
® JHE], we may choose now the smallest transition probabilities compatible with the 
current (fT3|). which are 

p /A , ,A x ^t\(P(q)~U*P(q)U)\^t) + h ,, 
MQ t+ r *q\Q t = q) = mP{q) \ %) • < 14 ) 

This need not coincide with the transition probability (jl2j) of (Qt), even though in the 
limit r — > also (Qt) converges to Bell's process. The same construction can be applied 
to the case \Q\ > 2 if U involves only pairs of configurations, i.e., if there is a partition 
of Q into subsets, all of which are either pairs or singlets, such that P(q)UP(q') = 
whenever q and q' do not belong to the same subset. Then (jlUJl still satisfies (JHJ) and thus 
defines a process. An example of this is a quantum computing circuit, realized through 
a time sequence of single qubit unitaries and CNOT gates. (Here, a configuration q 
corresponds to a definite value for the computational basis observable for each qubit.) 

To contrast the previous processes with an example of a process that does not con- 
verge to Bell's process in the limit r —>■ but has the quantum distribution (j3J) at every 
time, we define the process (Q1)terZ by the transition probability 

HQt+r = q\Qt = 4) = (% + r\P(q)\% +T ) ■ (is) 

This means that for every t, Ql is independent of the past and has the quantum distri- 
bution. Its limit as r — > 0, in a suitable sense, is simply the process (<3t)ieR f° r which 
every Q* t is independent of the past and has the quantum distribution, a process reminis- 
cent of Bell's j3] description of a precise version of the "many worlds" interpretation of 
quantum mechanics: "[Ijnstantaneous classical configurations [Q] are supposed to exist, 
and to be distributed [...] with probability \ip\ 2 - But no pairing of configurations at 
different times, as would be effected by the existence of trajectories, is supposed." 
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